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We follow recent formulations of dimensionally reduced loop operators for quantum field theories 
and exact representations of probabilistic lattice dynamics to identify a new scheme for the evaluation 
of partition function zeroes, allowing for the explicit analysis of quantum critical phenomena. This 
new approach gives partition function zeroes from a factored quantum loop operator basis and, 
as we show, constitutes an effective mapping of the renormalization group /3-function onto the 
noncommuting local operator basis of a countably finite Hilbert space. The Vafa-Witten theorem 
for CP-violation and related complex action problems of Euclidean Field theories are discussed, 
following recent treatments, and are shown to be natural consequences of the analyticity of the 
limiting distribution of these zeroes, and properties of vacuum regimes governed by a dominant 
quantum fluctuation in the vicinity of a renormalization group equation fixed point in the infrared. 
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A classical spin system and its quantum system ana- 
logue apparently share many similarities at the level of 
the spin Hamiltonian. The difference between quantum 
and classical formulations is perhaps somewhat subtle 
at first glance. In the latter case the spin components 
are represented through the action of an operator basis, 
rather than purely orthogonal vector basis elements as is 
the case with the former. Although the dimensionality of 
a particular model is unaffected by this step of operator 
inclusion the partition function of the quantum system is 
formally defined as a trace over this operator basis act- 
ing on some suitably structured Hilbert space (whereas 
the classical partition function is strictly given simply 
through a summation over spin vector realisations). This 
inclusion therefore adds an extra dimension to evaluation 
of the quantum partition function that maybe then dealt 
with, or expressed, through a variety of scenarios relating 
to the applicability of dimensional reduction. This basis 
reformulation or inclusion of operator components can in 
fact be argued to be a form of second quantisation of a 
given classical model Q. 

We are ultimately interested here to derive an exact 
polynomial expansion for the partition function for a 
generic quantum system formed in such an operator ba- 
sis, and to discuss the effect that quantum fluctuations 
have on frustrating this effort. The zeroes of this poly- 
nomial expansion then constitute the partition function 
zeroes whose limiting distribution we aim to now analyse 
following the arguments and reasoning presented in 
However, the difference from these two cases we will ar- 
gue (and would hope to stress as a central result) is that 
in general we identify the partition function of a generic 
quantum system is nonanalytic and unamenable to such 
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expression. Our reasoning will follow much as for the 
general case of renormalization group equations of mo- 
tion and the /3-function, which may become singular prior 
to reaching its asymptotic value through infrared diver- 
gences 0. Crucially dependent on structuring the for- 
mal limits in the countably finite Hilbert space basis we 
will argue a limiting analyticity can be defined, however, 
using recent results and discussions of the structure of 
the poles of the inverse Laplace transform in the con- 
text of recent results on exact lattice probabilistic dy- 
namics 

It may initially appear some what confused to attempt 
to define bounds for the analyticity of this expansion, or 
indeed of little practical relevance to order and arrange a 
form for an exact expansion that is valid essentially only 
in the vicinity of the continuum limit. This is certainly 
somewhat of an anathema to the scaling arguments of the 
renormalization group and partition function zeroes for- 
mulations where scale invariance is a central point. Our 
major thrust, however, will be to treat disconnected con- 
tributions to the partition function in a local basis, mark- 
ing a certain departure from (quasi-classical) perturba- 
tive expansion. The noncommutativity of the operators 
for the quantum system force the reformulation in this re- 
gard. Diagrammatically disconnected contributions will 
be here suppressed through the limiting distributions of 
the basis operators formed on the Hilbert space, rather 
than through the tacit convergence of any effective se- 
ries. Akin to partition function zeroes formulation, sin- 
gularities are expressed through the limiting asymptotic 
behaviour of poles as the continuum is approached, but 
here with an intermediate result closely connected with 
dimensional reduction that factorises a local symmetry 
preservation. To facilitate this discussion we first review 
salient points in the derivation of the transfer matrix and 
associated Quantum Monte Carlo method for assigning 
exact probabilistic dynamics to a loop operator basis, al- 
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though we would stress this is not primarily intended as 
a discussion to motivate specific numerical approaches. 



I. NONCOMMUTATIVITY OF THE 
HAMILTONIAN OPERATOR 

A number of algorithmic scenarios exist for the 
Monte Carlo simulation of quantum spin systems 
013 03 EH- They were developed for the purpose of 
addressing the critical slowing down that affects the sim- 
ulation of quantum systems at low temperatures. They 
differ largely from standard lattice gauge theory simula- 
tion methods, for example, in that in the main they are 
driven locally. Probabilistic Monte Carlo decisions are 
taken on a lattice site-by-site basis rather for the global 
lattice configuration as is usually the case with the lat- 
ter [l^. Although, we will later comment on the rele- 
vance of the results to the complex action problem and 
recent lattice attempts at a general solution in this di- 
rection here algorithms for the numerical study of quan- 
tum spin systems serve as a means to introduce the for- 
mulation a complete-orthonormal operator basis acting 
on a Hilbert space in the context of exact lattice prob- 
abilistic dynamics. Following recent discussion on the 
limiting continuum distribution of a suitably arranged 
Markov process, we then argue this is amenable to ana- 
lytic expression. Thus a means is identified to formulate 
partition function zeroes to express singular behavior in 
a quantum system for cases where otherwise the non- 
commutativity of the quantum operator will basis spoils 
this analyticity. This result is so expressed through an 
approach closely connected to dimensional reduction 0. 
From this derivation it becomes then possible to the de- 
fine bounds under which a generated Monte Carlo dis- 
tribution becomes suitably close to the continuum result 
that the zeroes of the partition function can be meaning- 
fully associated with the limiting thermodynamic argu- 
ments for the singularities of [2j |3j ■ Or more properly now 
their quantum system analogues, as we will argue from 
the results on the behaviour of the poles of the inverse 
Laplace transform of this expansion. In this sense, there- 
fore, determining the statistical limitations of numerical 
Monte Carlo sampled loop operator distributions with 
relevance to this task remains an onerous and somewhat 
separate ongoing direction, although we return to discus- 
sion of the complex action problem in a later section. 

Proceeding with commenting on generic Quantum 
Monte Carlo approaches for quantum spin systems, 
the pedagogical example usually taken is the spin one 
half anti- ferromagnetic Heiscnbcrg model (S = 1/2 
AFM)00. The transfer matrix of this system con- 
tains the simplest off-diagonal singlet contribution pos- 
sible to express a quantum interaction. The spins de- 
fined Si = 1/2 o \ are defined for the model on each 
lattice site index i of a d-dimcnsional hypercube, with a 
nearest-neighbour interaction and &i being the Pauli ma- 
trices. In general, ie. for all but vanishing off-diagonal 



contributions of the suitably defined transfer matrix, the 
spin operators do not commute. If one were to wish to 
evaluate the Boltzmann factor exactly from the nearest- 
neighbour spin Hamiltonian thus defined on a finite lat- 
tice volume this would then present something of a prob- 
lem, akin to the difficulties of representing Grassmann al- 
gebras in lattice gauge field theories. The combined idea 
of Suzuki and Trotter as a basis for numerical simulation 
was to partition the nearest-neighbour Hamiltonian op- 
erator H = Y.i H i- For tnc case 01 the 5=1/2 AFM the 
leading result would be then given for an odd-even site 
partitioning of the operators. 

e-rf «[]exp(-Ar J ff l ) + C , (AT 2 ) (1) 

i 

The error in estimating the exponential is thus sub- 
dominant despite the noncommutativity of the nearest- 
neighbour spin operators that constitute the Hamiltonian 
7i. This is dependant evidently on the value of At, as 
yet undefined. The partition function itself is defined as 
the trace over the basis on which the operators act, and 
thus the extent of the above product i must clearly be 
somehow dependent on the lattice volume of the classi- 
cal spin model on which the quantum system operators 
are introduced. For the 5=1/2 AFM, as we have said, 
the leading result is given through odd-even partitioning. 
There would be no reason a priori, though, to define a 
particular lattice extent for the additional dimension in- 
troduced by the operators themselves. Indeed, only two 
points need be formally defined for the trace operation. 
If At = (3/L t is the defined as the lattice spacing in Eu- 
clidean time the only necessary points to preserve the 
measure are At = and At = j3. However, we wish 
to discretise the Euclidean time direction, not simply for 
the purposes of numerical simulation but to now enable 
discussion of the recent limiting results of probabilistic 
dynamics on this Euclidean extent by first treating it as 
a countably finite space. To properly define the partition 
function Z(/3), it is necessary to first note that the ele- 
ments of the factorised product are the transfer matrices 
that act on the suitably arranged complete-orthonormal 
set of basis states. 

Z((3) ee Tre"^« ^n«) ( 2 ) 
{S u } P 

Wp(Sp) = {Si t tSi+i t t\e T ~ H *|iS, j t-j_i 1 S,_|_i i t-|_i) (3) 

The summation of the weights W p for the transfer matri- 
ces above therefore runs over the spin configuration space 
of the lattice, as given through the state of the elemen- 
tary plaquettes p. For the 5=1/2 AFM these are defined 
as the two-by-two elements of the hypercube defined over 
a unit spatial i and Euclidean time t extent. The weight 
W p associated with the of each plaquctte are given by the 
matrix elements, and can be built up numerically through 
the probability density function distributions generated 
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by a Monte Carlo procedure, as we will later discuss [Icj . 
The definition of the transfer matrix elements and their 
explicit relation to the noncommutativity of the spin op- 
erators will also be postponed to later section, whilst we 
consider an alternate quasi-classical approach [llj] 0| ■ 

The reverse of the above product expansion is of 
course possible. Rather than proceeding from a lo- 
cal, and inherently limiting-case factorisation perspec- 
tive, one can aim to treat the statistical fluctuations of 
the noncommuting operator basis in a global manner, us- 
ing the ideas of stochastic field theory. Statistical cumu- 
lants (Hi . . . Hn) c can be defined on the operator basis 
through the differentiation of some suitable generating 
function, 
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It is important to note that no assumption is being made 
for the distribution of the operators on the space in this. 
If this were the case one would arrive at a standard time- 
ordered path integral formalism for quantum field theory 
generating functionals, albeit on a countably finite space 
(a Fock space definition), rather than this statistical re- 
sult which is simply for functions rather than functionals. 
That the time-ordered result is not guaranteed is some- 
how obvious since the operators in the cumulant defini- 
tion are as yet rather general and undefined. Their action 
could be highly nonlocal or could indeed be incomplete 
in some sense, in which case the dynamics would be beset 
by singularities. That the action of a proper Hamiltonian 
necessarily yields a time-ordered product is somehow less 
obvious, but will we again postpone a fuller discussion of 
issues of locality and noncommutativity to the section 
that follows, where the role of a Markov process in defin- 
ing exact probabilistic dynamics for the lattice is made 
evident. 

Following this analogy, though, the expectation value 
of the product defined without the logarithm in the right- 
hand side would correspond to analogue of the mixture 
of disconnected and connected n-point functions one ob- 
tains from the action of the generating functional in the 
field theory analogue. Again, without assuming form for 
the distributions one can take the continuum limit of the 
operator space N — > oo for the above definition of cumu- 
lants. This then yields a quasi-classical effective Hamilto- 
nian Ti. c from equating the cumulant definition to an ex- 
act series expansion for an exponential in this limit. This 
appears evidently simpler to evaluate, being reduced to 
a quasi-classical formulation and therefore amenable to 
a summation over vector realisations. 



Z= /n(n|exp(-/3H)|n) (5) 



Caution is in order however as one must now expand in 
the matrix elements of powers of the full Hamiltonian H, 
for a practical numerical approach to evaluate the quasi- 
classical Hamiltonian TL C . 



0H c (n,[3)= {n\H\ny-^(n\H 2 \n) c 



3! 



(n|^|n) c + ... (7) 



exp[-0Hc(n,0)] = (n\exp(-/3H)\n) (6) 



As was case with the Trotter-Suzuki approach the non- 
commutativity of the basis cannot be entirely circum- 
vented. This expansion may appear somewhat straight 
forward, but H c itself is strictly defined only through the 
cumulants (connected components), rather than the sta- 
tistical expectations (connected and disconnected com- 
ponents). A numerical evaluation of the quasi-classical 
system following similar truncation approximations to 
the Trotter-Suzuki form therefore suffers from three 
sources of systematic error. The expansion itself may di- 
verge through either the value of (3 or the level-spacing di- 
verging as a function of the expansion, and the stochastic 
expansion terms themselves contain the statistical ana- 
logue of disconnected contributions to the quasi-classical 
partition function in all but the continuum limit. With- 
out unduly pessimistic focus on the numerical prospects, 
both of the two parallel approaches fall on this issue of 
noncommutativity. Yet at the same time, this compari- 
son as the continuum limits are separately approached is 
revealing. 

To make a further association with issues of dimen- 
sional reduction in gauge field theories we can aim to 
relate quasi-classical statistical and functional terminolo- 
gies for effective actions. At small inverse temperature 
P, the interaction coupling of gauge fields becomes small 
through asymptotic freedom and perturbative expansion 
is justified For length scales much larger than 

those of the Euclidean time extent it also seems rea- 
sonable to integrate out the nonstatic (time-dependent) 
modes of the gauge fields. This gives the leading behav- 
ior at small f3, but it is important that the perturbative 
expansion of this dimensionally reduced model does not 
have divergences in the infrared. The nonstatic modes 
are suppressed and renormalize the expansion terms of 
the new effective theory. To develop the above statistical 
association, the expansion is again valid if the perturba- 
tive coupling is not so large as to cause divergence. By 
way of contrast no explicit assumption is made about 
the locality of the expansion in the statistical case and 
so it is not renormalized. The dimensional reduction of 
the field theory here forms the functional analogue of the 
stochastic approximation. The functional approach be- 
comes again invalid if the contribution of higher order 
n-point functions (cither connected or disconnected) be- 
comes dominant. 

It is essentially the locality of the expansion of the ef- 
fective action that crucially determines the response to 
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infrared divergences and so we seek a different approach 
as a basis for study. By following the Trotter-Suzuki de- 
composition, as opposed to the effective quasi-classical 
Hamiltonian formalism, we work in a nonlocal frame- 
work without implicit effective series truncation. Lee 
and Yang, in the context of classical spin systems de- 
fined a limiting analyticity for the partition function such 
that thermodynamic singularities were only strictly ex- 
pressed through the limiting distributions of the zeroes 
of polynomials so constructed ■ We we will argue that 
this specific analyticity property of the partition function 
is a nontrivial consequence of the renormalization group 
approach in the quantum case by working over this ex- 
tended phase space. By addressing the form of the limit- 
ing distribution taken by noncommuting operators on the 
basis space we will seek to similarly define a limiting an- 
alyticity for the infrared divergences connected with the 
renormalization group equation fixed points of quantum 
systems through this factorisation approach. 



II. LOOP-CLUSTER OPERATORS 

Before proceeding to discussion of the partition func- 
tion polynomial form for the quantum case we con- 
tinue first outlining the relationship between the fac- 
torised Hamiltonian operators and their distributions on 
a Hilbcrt space. This is the relationship defined through 
the transfer matrix elements of these operators. We moti- 
vate these now necessary definitions through part-review 
of numerical Quantum Monte Carlo loop-cluster meth- 
ods. Numerical motivation perhaps seems an unneces- 
sarily heuristic path to take here for the proposing of an 
exact continuum form for the generic partition function 
of a quantum system. To reiterate, indeed the stated aim 
is not to define numerical grounds for simulation. By con- 
structing, or reviewing, the Markov process definition we 
will simply articulate the complctc-orthonormal nature 
of the basis formed by the elementary plaqucttes defined 
over the spatial and temporal sites, and exact probabilis- 
tic dynamics defined thereon. 

Again whilst a factorisation property is clearly a key 
clement in properly defining the zeroes of a polynomial, 
the direct connection with a complete-orthonormal ba- 
sis for a Hilbert space and partition function zeroes of 
a generic quantum system would perhaps appear obtuse 
at this point. However, mirroring Lee and Yang's dis- 
cussion we associate the analyticity of an exact expan- 
sion with identifying the continuous distribution or func- 
tional assignment of solutions to the generating functions 
of said expansion. That each expansion coefficient in the 
polynomial expansion can represent a well-defined inte- 
gral equation is actually a necessary starting point of Lec 
and Yang's discussion Q ■ We therefore do this explicitly 
here by considering the action of taking the vanishing 
Euclidean time lattice spacing limit for the loop opera- 
tors on a finite spatial lattice extent. The construction of 
the quantum system analogue of Lee and Yang's integral 



assignment is then given through this action, implicit to 
the loop-cluster method definition. 

The loop-cluster method for quantum systems is simi- 
lar in form to the that of the Swendsen-Wang cluster al- 
gorithm of classical statistical mechanics |l7j [7] 9] . The 
partition function is defined over a set of spin configu- 
rations {S} of the lattice, and a set of graphs {G} link- 
ing these spin states, defining an extended phase space. 
For the Monte Carlo process statistical weight functions 
W(S, G) must be assigned for each spin configuration of 
a lattice ensemble, such that now all elements of the ex- 
tended phase space are defined as cither being connected 
or disconnected. Given some configuration of the spins 
S of a lattice each bond linking spins is separately con- 
sidered and assigned such a status based on the value of 
W(S,G), determined through the local spin configura- 
tion of S. This then implicitly defines a path, or graph 
subset G, through the volume of the particular configu- 
ration. The total weight of a given configuration W(S) 
in the lattice ensemble is given by the sum over graph 
weights. 



J2W(S,G) = W(S) p(S^(S,G)) = -^y, 

G 

(8) 

where p( S — ► (S, G) ) is probability of a obtaining a par- 
ticular graph within a lattice configuration. This then 
implies both that W{S) ± and W(S, G) £ 0. The 
process is otherwise noncrgodic and there is some abso- 
lute preferred configuration state of the volume that once 
reached cannot be left in an exact statistical sense. We 
are interested to treat singularities, of course, but do not 
anticipate divergences such as these on a finite lattice for 
essentially geometric reasons - without, though, wishing 
to necessarily restrict the spaces of interest at present. 
The loop-cluster algorithm itself is clearly of more gen- 
eral construction than its applicability to quantum sys- 
tems, but the Hamiltonian operation we are interested 
in, acting locally on suitably defined plaquettes, can sat- 
isfy this latter Monte Carlo probabilities detailed balance 
condition. 

Returning to the importance of locality. Considering 
the pedagogical example of the 5 = 1/2 AFM according 
to the above definition of a graph a single spin cannot 
be inverted alone on an elementary plaquette to create 
a valid change between successive lattice configurations. 
This follows from the implicit local conservation of S. 
The graphs so defined thus map relative change across 
the extended volume. In fact, demanding on a given a 
graph that the form of W plaq (S, G) does not change at all 
between consecutive lattice configurations in an ensemble 
then implies that all allowed plaquette changes match up 
between neighbouring plaquettes across the volume. 

An extensive literature exists on the definition of suit- 
able definition of transfer matrices for quantum systems 
and the loop-cluster method, but to clarify further for the 
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pedagogical example the here plaquette has eight spin 
configurations, S p = i^, of the four site vertices, i Q- 
Each transition i «-> j defines a respective transfer ma- 
trix element w lJ = w Jl , and corresponds in this case to 
inverting either two or all four of the spins of a plaquette. 
We can then therefore think to write the partition func- 
tion for this, or a generic quantum system, in terms of an 
equivalent bond or loop operator basis defined onl y in 
terms of the pertinent transfer matrix elements Q 116^|18|| 
. The pertinent transfer matrix elements will be those are 
nonzero. The diagonal elements of these loop operators 
for the pedagogical case will correspond to four inversions 
and comprise the identity operator, and the self-adjoint 
off-diagonal remainder will be given uniquely through the 
graphs determined for each configuration (or rather their 
ensemble-average). This all then essentially represents 
an effort to define the pseudo-reduction step on a finite 
basis, mapping the discrete distribution of graphs to the 
limit of vanishing Euclidean time lattice spacing where a 
continuous distribution is recovered. The partition func- 
tion in terms of the loop operators h b is, 

Tre"^ = e ^^ A lim o (jJ b e(- J+J ^ A *)^ /A ' (9) 

= e?V*J j p(duj)f[h b , (10) 

where the Hamiltonian is rewritten over the bonds b as 
H = — J2b J hb and J is the nearest neighbour coupling. 
Here the Y[* hf, is a time ordered product of the loop oper- 
ators with p(dio) a Poissonian probability measure of the 
connection status recovered by taking the limit of vanish- 
ing Euclidean time lattice spacing. That surprisingly this 
relatively general construction of the loop operator trans- 
fer matrices yields a precise continuous distribution in the 
limiting step follows directly from the locality property. 
It is important to note as well that no noncommutativity 
property is invoked to obtain the Poissonian distribution 
on the basis elements of the now countably finite spatial 
element space. One could have equally well ignore the 
off-diagonal elements altogether and used the approach 
to map a classical system into another classical system, 
since after all the method is classical in origin. 

The advantage over the quasi-classical series expansion 
approach here is that the reduced factorised basis ele- 
ments are explicitly highly nonlocal in character whilst 
locally preserving S exactly. The central assumption, 
however, is again that of the quasi-classical expansion 
: the statistical averages only become interchangeable 
with those of the thermodynamic or dynamic averages in 
the infinite volume limit. Here the definition is somehow 
less ad hoc and the treatment of disconnected contri- 
butions less tacit, since locality is enforced site-by-site, 
which allows us define a polynomial with a slightly dif- 
ferent relation to the poles of the partition function. By 
construction, the basis of plaqucttes is complete, and or- 
thonormal and more importantly the weight functions 
associated to graphs are unique and independent of the 



location within the Markov chain or statistical ensemble. 
Although S is locally preserved the elements of the ex- 
tended space reached by a finite Markov process can be 
necessarily incomplete and so functional solution is not 
necessarily admitted. We could be trite and label this the 
ambiguities of numerical determination, but the issue of 
whether or not a discrete extended basis can form a cov- 
ering of the target space is a perhaps a more interesting 
question which we will return to shortly. 

The extent to which the above partition function corre- 
sponds with the Nambu definition of a generalised Hamil- 
toni an ( in the operation of elements on a Poisson mani- 
fold [n|), and could perhaps then therefore be now fur- 
ther generalised for non-Euclidean metrics, is certainly an 
interesting question but not one we will aim to further 
address here. We wish to now investigate the limiting 
analyticity of this factorised reduced expansion, noticing 
both that the probability distributions of the operators 
become exact and continuous in the limit of vanishing 
Euclidean time lattice spacing and that a product of two- 
by-two self-adjoint operators (for the S = 1/2 AFM case, 
and of generic form) is particularly simple to diagonalise. 
We have now identified the continuous distributions as- 
sociated with factored product which we are aiming to 
now diagonalise, in keeping with Lee and Yang's motiva- 
tions, as a first step towards establishing the analyticity 
properties of the eigenvalues and characteristic polyno- 
mial associated with this operator product. 



III. PARTITION FUNCTION ZEROES 
POLYNOMIAL 

We have stressed, perhaps seemingly unnecessarily, 
that the partition function zeroes for a generic quantum 
system are somehow different than those introduced by 
Lee, Yang and Fisher which we interpret now as classi- 
cal results. By considering the action of loop operator 
elements on an extended phase space, defined through a 
complcte-orthonormal local basis, we have discussed in 
the previous section the identification of an exact fac- 
torised operator product on a countably finite Hilbcrt 
space with an associated continuous time-ordered distri- 
bution structuring. This is a limiting result, but as we 
have noted the finite basis loop operators by construc- 
tion are trivially diagonalisablc in general. The new step 
we now take is to note that there is a necessarily a differ- 
ence between diagonalising the elements of a direct prod- 
uct and determining the characteristic equation of that 
product. Herein, we believe, lies the divide in the argu- 
ments concerning deriving an exact analytic form for the 
partition functions of the respective quantum and classi- 
cal systems. 

If the distribution obtained through projecting the 
graph results onto the spatial extent of the lattice is dis- 
crete, rather than continuous, in the quantum case some 
of the loop operator elements may be defined only on 
a subspace of the extended space. The transfer matrix 
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allows for the lattice to be entirely comprised of discon- 
nected contributions, of course, but this singularity is un- 
der defined in the discrete quantum case. The eigenvalues 
of the factored product are therefore not necessarily an- 
alytically connected on the finite Hilbert space defined 
over the spatial elements with a discrete statistical dis- 
tribution, ie. at finite Euclidean time lattice spacing. 
Conversely with what we here deem the classical result 
(if this pictured simply as a limiting case of commuting 
operators) all elements are necessarily connected through 
the diagonality of all the terms of the factored product. 
There may be still disconnected contributions to the par- 
tition function but the space of the partition function 
itself is not now under defined as a consequence, since 
the plaquettes are constructed by definition to preserve 
a local symmetry Up to singularities in the continuum 
limit for the Hilbert space, as defined over the loop op- 
erator basis of the lattice spatial extent, we can thus be 
assured that the classical case is always defined whereas 
in the general quantum case is not. At least, to be clear, 
if we arc now aiming to treat the continuum limit of the 
eigenvalues of our product as the polynomial zeroes of 
the partition function of Lee, Yang and Fisher we must 
treat the spatial and Euclidean time extents separately 
in our analyticity arguments. This follows essentially as 
a direct consequence of the noncommutativity of the op- 
erator elements of the Hamiltonian in the quantum case. 

In the quasi-classical series expansion formalism we 
would now have difficulty proceeding with our discussion 
as a basis for practical calculations because the basis op- 
erations in the extended phase are nonlocal. This means 
that we could only recover both of our limits simulta- 
neously : the discrete statistical distribution being built 
up step-by-step is that for the entire volume, rather than 
for the elementary plaquettes alone. Conversely, the fac- 
tored product elements we will now focus on are being 
defined towards the limit of vanishing Euclidean time lat- 
tice spacing, thus giving an analyticity ahead of the con- 
tinuum limit for the spatial extent through a preserved 
local symmetry. It is this effective treatment of the dis- 
connected contributions that allows us to proceed in the 
quantum case with the factorised product approach. 

What we want to now consider, therefore, is if the 
eigenvalue spectra, assumed now complete and cover- 
ing for the extended space, has a meaningful connection 
with the poles of the free energy density of the parti- 
tion function and if these poles then demonstrate the 
lattice volume scaling properties associated with parti- 
tion function zeroes analysis of Lee, Yang and Fisher. 
We do this now by noting recent results on the struc- 
ture of the inverse Laplace transform poles of the trans- 
fer matrix elements in the context of discussions on ex- 
act probabilistic lattice dynamics. To follow this recent 
discussion, and deduce our zeroes continuum limit be- 
haviour, we can rewrite the loop operator expression de- 
termined for the partition function through the notation 
of 010 I n this case multiplying by the infinitesimal 
Poissonian probability elements in the Euclidean time 



intervals [ti, ti + dti], . . ■ , [tjf, ijv + dijv] the off-diagonal 
operator matrix elements WV (/?) of hb are now expressed 
as, 

Wn(P;V ,Vi,...,V n ) = 

e N jjj £ dU 0(U - cxpAA)| e- v ^{ll) 

with Aj = Vi — Vi-i the potential energy difference be- 
tween the sites at i and i— 1, and e some arbitrary energy 
scale for the purposes of normalisation. In the previous 
section we showed that in the limit of vanishing Euclidean 
time lattice spacing we recover a time-ordered Poisson 
distribution for the elements of the factored spatial basis. 
Here this is the initial assumption for the product : an 
element remains in a given basis state until it undergoes a 
Poissonian distributed event. We can thus consider this 
the generic form for the loop operator elements in this 
limit. Here we are again expanding the product over the 
basis of states N, the extent of which corresponds to the 
spatial volume, with the integrand here corresponding to 
the time-ordered limit of the graph distributions. In the 
discrete case the graph paths can propagate throughout 
the configuration volume, and equally here the product 
elements are in causal connection. The above element 
maps between consecutive spatial sites i, with all four 
corners of the elementary plaqucttc considered by defini- 
tion, and again S is locally conserved. The explicit dif- 
ference, however, between this weight form Wjv(/3) and 
hb lies in the expression of the specific matrix element 
contributions. The Heaviside step function being uni- 
directional over the product indexing here implies that 
the mapping of Wjv(/3) defines only one element of the 
transfer matrix. The time-ordering is here implicit to 
the potential energy values recovered from the graph de- 
termination process, whereas the ordering of the product 
indexing simply labels the element within the elementary 
plaquettc. Since this latter ordering maps between suc- 
cessive site indices i we can associate Wn (P) with the off- 
diagonal contribution to the loop operator hb- Why we 
believe it relevant to consider just this subset of hb is that 
for the generic form of hb for a quantum system, as we 
have argued (somewhat tacitly) in the previous section, 
the diagonal elements form the identity operator. The 
zeroes of the operator product are thus characterised by 
the values of the self-adjoint off-diagonal contributions. 

Notice again that the above expression is a result for 
an exactly ergodic Markov process, and so we must have 
achieved this statistical limit of vanishing Euclidean time 
lattice spacing for the above relation to be valid. We have 
argued, somewhat hcuristically, that this limit is reached 
by generating a sufficient number of the formally defined 
graphs such that these graphs then recover a suitably 
complete basis for the extended space. Since the trans- 
fer matrix elements are as well defined over the ensemble 
rather than for a single lattice configuration, similarly we 
will be interested in the zeroes of the ensemble-averaged 
characteristic polynomial rather than those of a single 
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configuration for a finite Euclidean time extent. The 
above relation for the weight form Wjv(/?) is of course 
given in terms of potential energies rather than the poly- 
nomial zeroes of the loop operator product, but the an- 
alyticity of the zeroes and poles are in fact intimately 
connected through the properties of the inverse Laplace 
transform of Wi\rQ3), which we will now consider follow- 
ing [j| [|| . Differentiating with respect to (3 and substitut- 
ing for the Laplace transform Wat(z) of Wn(P) reduces 
the weight form to the following recursive equation, 

/>oo 

Wn{P)= dl3e- zl3 W N (P), z&C (12) 
Jo 

dpW N (f3) = eW N -x(0)-V N W N {P), (13) 
W N (z) = e{z + V N )- x W N -i{z) : (14) 

- e ~n^. do 

i— 

The motivation for introducing these new relations is es- 
sentially to quantify the asymptotic scaling behaviour of 
the expansion via discussion of saddle point solutions. 
The above poles of the Laplace transform are by defini- 
tion confined to be in the real segment [— V m ax, — Knin] , 
with Vi nax and V m in respectively the maximum and min- 
imum elements of the set of potential energies {Vi}. It 
follows from the definition of the Laplace inverse trans- 
formation that its integration contour T is any line par- 
allel to the imaginary axis contained in the analyticity 
domain of the Laplace transform, ie. Rez > — Vm m . Ap- 
plying a saddle point approach to the expansion it can 
be shown that the first N — 1 solutions lie in the range 
— V max < Rezo < — Vmin with the last solution lying such 
that Rezo > -Vmin , where z are the locations of the 
saddle points and the ordering of solutions is in terms of 
increasing size. The restrictions on the poles come from 
quite general considerations, the only real condition here 
being that j3 is necessarily positive [23 [2l| ■ 

The whole point or purpose of the saddle point anal- 
ysis here is that its solution is independent of the path 
distributions of the graphs but rather dependent on the 
local properties of the graphs as given through the po- 
tential energy differences. From this consideration then 
the integration contour can be deformed to a new par- 
allel one, although we must know by definition that the 
path mapped by the weight is analytic along the con- 
tour r and where the isolated singularities exist. One 
could therefore then in principle equally well treat the 
general case of complex-valued potential energy differ- 
ences, which would be similarly defined since L would 
have no explicit singularities from the local properties of 
the graphs by construction. Generally, of course, there is 
no reason to suppose that the poles have to be uniquely 
defined, and could be given only up to a branch cut, but 
this is not a problem with the asymptotic form since in 
the asymptotic limit this choice is naturally resolved. We 



should note that through the above definitions and those 
of the previous section we are at the same time uniquely 
establishing the residue properties of the weight function 
and thus the existence of the zeroes of the loop operator 
factored polynomial. 

The last saddle point solution alone is compatible 
with the restrictions on the integration contour T, and 
can be further shown to be bounded from below by 
Rez > -Vmi„ + N Vml J/3 0. Thus for (3 -> oo with 
Nv mln ~ N ~ (3 the distance between the saddle point 
and the closest pole remains finite, and the zero temper- 
ature continuum limit is well-defined. The distribution 
of the properly defined complete zeroes now necessar- 
ily have some limiting distribution with relevance to the 
continuum from the asymptotic behaviour of these poles. 
To be clear, the partition function is expanded in what- 
ever constant diagonal prefactor is introduced in hb- It is 
proportional to /3 for the pedagogical S =1/2 AFM ex- 
ample and therefore the analysis is specifically of interest 
for (although not restricted to) the case of quantum fluc- 
tuations at low temperatures. Strictly, as we have seen, 
only one of the poles lies in the analyticity domain of the 
polynomial. This does not imply that the characteristic 
polynomial is undefined, but rather that we only under- 
stand the asymptotic behaviour of the weight function in 
the vicinity of this one pole. There is a circle theorem 
by Lee and Yang which we will now discuss which relates 
to the limiting behaviour of branch cuts, which we can 
perhaps envisage now being applicable in at least one in- 
stance of the zero temperature limit. To reiterate, the 
noncommuting nature of the quantum system Hamilto- 
nian has forced us to consider a limiting Euclidean time 
spacing distribution defined over the spatial basis of our 
lattice, for which at least one of the zeroes of the polyno- 
mial by construction has a well defined continuum limit. 



IV. SCALING OF PARTITION FUNCTION 
ZEROES 

We have considered how a partition function can be 
defined over an extended phase space for a quantum sys- 
tem through a loop operator formalism, and also how 
the limit of vanishing Euclidean time lattice spacing is 
related to the analyticity properties of the characteristic 
polynomial of this factored product and its pole struc- 
ture. Although the recurrence relations defined for the 
above Laplace transforms of the operator elements have 
an explicit dependence on the spatial lattice volume N, 
we must further explore this connection to firmly estab- 
lish the zeroes content regarding the singularities of the 
partition function as the limit N — > oo is approached. The 
starting point of Lee and Yang's arguments is essentially 
that a valid zeroes polynomial form exists, from which 
finite volume scaling discussion can be then made f|. 
We have now argued that the analyticity of factorisation 
property of a quantum partition function is in general not 
possible without first approaching an asymptotic scaling 
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regime : a far weaker precondition for this zeroes form. 
Comparing with the constraints we have thus far defined 
on the factored product eigenvalues for the quantum sys- 
tem we now contrast with the two central criterion of Lee 
and Yang's discussion for defining valid thermodynamic 
functions from the partition function in the limit ./V — > oo 
of the spatial lattice volume 0] ■ 

The first of these criterion is that the limiting value of 
as the pressure of the system P, is independent of the 
shape of N and is a continuous monotonically increasing 
function of the expansion parameter of the suitably fac- 
tored partition function. 



P = lim 



1 

N 



log(Z) 



(16) 



Since the action of the factorisation in the quantum 
case is by definition local and ergodic it follows that the 
weights are obtained recursively, as we have seen. The 
limit (TV + M) ^ oo of the difference of the log of two con- 
secutive intensive weights of integer index N and (N+M) 
is thus zero. Establishing this property one can then con- 
sider the volume as consisting of a finite subset of smaller 
volume elements, with interactions across the boundaries 
between subsets for which the free energy of this inter- 
face tension does not grow faster than the surface area. 
The above result thus implies that the limit of thermody- 
namic functions, including the pressure, is always finite. 
Monotonicity and the continuous nature of the functions 
follows from the similar behaviour of the derivatives of 
the partition function, as is seen explicitly through the 
properties of the inverse Laplace transform. Thus the dis- 
tribution of the zeroes necessarily gives all the analytic 
behaviour of thermodynamic functions in the expansion 
parameter plane as N— >oo. 

The second of Lee and Yang's criterion is that the lim- 
iting value of the density of the system p, is an increasing 
function of the expansion parameter of the suitably fac- 
tored partition function but does not necessarily take a 
defined limit for all expansion parameter values, (3. 



lim 



d 



1 



n^oo dlog(P) N 



(17) 



Again this is very much an implicit result for the factored 
product of the quantum system Laplace transform and 
the discussions of the previous section . The criterion 
implies that if a region of the positive real axis of j3 is 
free of zeroes of the partition function as N — > oo then 
the system is here in a single phase. The proof of this 
follows from noticing that the coefficients of the factored 
product expansion are necessarily positive, which is a 
fairly general consequence of the operators being self- 
adjoint in the quantum case. 

Lee and Yangs first discussion with their new result was 
essentially to notice that there was necessarily a problem 
in discussing thermodynamic functions with solely the 



TV — > oo limit of the partition function polynomial expan- 
sion coefficients : this leads to the result of undefined 
phases. It was not possible to discuss thermodynamic 
singularities directly in the asymptotic limit, rather the 
convergence to that limit was the important consider- 
ation. They reviewed an earlier treatment of Mayer 
in which they showed it becomes impossible to analyt- 
ically continue the series expansion thus formed between 
phases. This is unlike of course, by construction, the 
above implicit spatial lattice volume dependencies of P 
and p given from the properties of the Laplace trans- 
forms of the operator elements that do yield well defined 
thermodynamic functions. The difference lies in the dif- 
fering properties of a series and polynomial. The series 
expansion values do not yield a density dependence to 
pressure for specific volumes below the point at which 
the zeroes approach the real (3 axis (the phase transition 
point). This is implicit since the continuum values have 
obviously no spatial lattice volume dependence. One can 
therefore argue from this that not only is the partition 
function zeroes form for the quantum system compati- 
ble with Lee and Yang's discussion, but also that it is 
therefore in some sense a necessary precondition for all 
partition function zeroes expansions, classical and quan- 
tum. To expand : we have laboured to ensure the fac- 
tored product is meromorphic, singular only at the poles, 
which arc themselves in a one-to-one correspondence with 
the zeroes. A loop operator picture defined over an ex- 
tended phase space must be therefore factorisable for all 
the thermodynamic functions to be able to defined for a 
system, since it is only the existence of the limit of van- 
ishing Euclidean time lattice spacing that ensures that 
the polynomial expansion zeroes, and by extension the 
expansion coefficients, have a valid asymptotic form. 

We can generalise this discussion further to consider 
the implications for complex operator elements and the 
introduction of symmetry breaking source terms defined 
with relation to the limit N — >oo. Partition function ze- 
roes have been recently used to discuss the Vafa-Wittcn 
theorem and complex action problems for quantum sys- 
tems [2^. To summarise this discussion : the Vafa- 
Witten result is that parity and CT cannot be spon- 
taneously broken in vector-like theories such as QCD. 
The argument follows from two points ; any arbitrary lo- 
cal hermitian order parameter for parity X, constructed 
from Bosc fields is proportional to an odd power of the 
totally antisymmetric tensor, and the free energy density 
is well-defined in the presence of a symmetry breaking 
source 8X, with 9 some real constant. These two con- 
ditions jointly then imply that the symmetry breaking 
source term be a phase in the integrand of the Euclidean 
partition function. The explicit connection of the theo- 
rem with partition function zeroes is then made by con- 
sidering p(X, N) the probability density function of X, 
where X expresses X extensively. Assuming parity to 
be spontaneously broken p(X, N) can then be argued to 
be taken as the sum of two delta functions in the limit 
N — >oo, centred on ±a. The relevant zeroes of the parti- 
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tion function in are thus obtained as the solutions of, 

C dX p(X,N) sm(8N(X - a)) 

cot(6Wa) = ) (- . (18) 

/ °° dX p(X, N) cosf 9N(X - a) 

This has an infinite number of solutions approaching the 
origin, and so the Vafa-Wittcn theorem is thus seen to 
be upheld since the second tenet of free energy is then 
undefined. In fact this is a general consequence of the 
properties of the Laplace transforms of complex-valued 
integrands. The asymptotic convergence, as a function 
of 9, is in general not defined for complex- valued func- 
tions given for the exponent : the oscillatory phase means 
the integration contour T is ill-defined. Conversely, as 
we have argued, following Lee and Yang's discussion of 
Mayer's scries treatment it is possible to consider the 
symmetry breaking consequences of a complex-valued 
system if the expansion coefficients are expressed inten- 
sively. Namely, if we have a contribution to the partition 
function that is highly nonlocal but preserves an exact lo- 
cal symmetry the asymptotic convergence properties of 
the partition function are defined, as we have seen ex- 
plicitly with the above loop operator formalism. 

We can understand in a very simple way how singu- 
larities of the partition function arise as N — > oo from 
our pedagogical loop operator example. Consider a self 
adjoint matrix with real elements of values between zero 
and one and the product of all elements in the upper tri- 
angular region such that these give recursively the char- 
acteristic polynomial expansion coefficients to successive 
order [23j. It can be shown by induction that if these 
elements have this probabilistic interpretation then nec- 
essarily the zeroes of the characteristic polynomial are 
bounded on the unit circle. In the case of Lee and Yang's 
discussion for the Ising model with external field we can 
then identify the field associated with generating the off 
diagonal transfer matrix elements and a symmetry break- 
ing nonzero source value. Although here there is no such 
association we can still understand the singularities of 
the spatial extent of the partition function in the same 
way, and therefore can similarly associate a critical value 
to our factored product expansion parameter with a valid 
asymptotic form. Generalising for complex valued oper- 
ator elements we see that the critical value of coupling 
is still recovered through a similar boundcdness induc- 
tion argument as the operators are necessarily hermitian 
in this case. We have in the loop operator introduction 
discussion and subsequent factorised product discussion 
essentially defined a source term such that it is only re- 
covered in the asymptotic limit, thus obviating the ambi- 
guity of contour definition for the Laplace transform for 
the complex case. 

Our initial motivations included a comparison of quasi- 
classical and factored product expansions for the treat- 
ment of quantum systems. It is perhaps now more ev- 
ident why we have favoured the latter for this zeroes 



treatment : the complex action problem above being 
amenable to treatment via quantum fluctuations. The 
equivalence of the quasi-classical approach to the con- 
tinuum is only resolved when both spatial and time ex- 
tent limits are simultaneously reached, but more over 
the symmetry breaking mechanisms are otherwise po- 
tentially contaminated with disconnected contributions 
to the partition function. Conversely, the nonlocal na- 
ture of the intermediate result of Euclidean time lattice 
spacing for the factored product case means that nonlo- 
cal symmetry breaking terms can be considered but on 
a countably finite space, yielding well-defined thermody- 
namic functions in the continuum. The quantum fluc- 
tuations represented through the partition function are 
inherently nonlocal and it is only through this latter lim- 
iting formalism that the disconnected contributions are 
exactly articulated since otherwise, as we have seen, the 
generating functional for the partition function becomes 
simply a generating function. 

There is an intrinsic connection between this discussion 
and the renormalization group arguments made by Wil- 
son through the implicit differentiability conditions intro- 
duced for the spatial lattice volume extent [24| ■ Looking 
at the density relation, the condition used to establish the 
N — > oo zeroes singularities, it is possible to understand 
that the critical points on the real axis of the expansion 
parameter correspond to the fixed points of the renormal- 
ization group equation defined in the limit N ~kx>. This 
is a new result of the loop operator treatment since for 
the general classical system partition function zeroes case 
it is not possible to make this association, as the asymp- 
totic scaling properties of the expansion parameter are 
not defined. The expansion parameter in this case can be 
strongly dependent on the Euclidean time lattice spacing 
and singular and dependent on the asymptotic properties 
therein. The differentiability conditions of the renormal- 
ization group equations have an obvious expression in 
the recursion relation for the Laplace transforms we have 
discussed, and the group property can be understood in 
this mapping between successive poles. Returning to the 
Vafa-Witten comment it is of course possible as with the 
renormalization group discussion that discontinuities be- 
tween minima of the vacua arise in the asymptotic limit, 
in which case it is possible to discuss phase transitions 
even in the case that we have no intensive source term 
comparing with the Ising model with external field |22j . 
The difference in level between the peaks of the vacua es- 
sentially then does not have to be infinite in order that an 
exact discontinuity arises in the asymptotic limit of the 
expansion parameter f3 which would now play an analo- 
gous role to time in the motivation of the renormalization 
group equations. 

In this loop operator treatment it is therefore impor- 
tant to understand the evolution of the spatial lattice 
spacing dependence in regards to the nonlocal contribu- 
tions to the expansion : infrared divergences potentially 
lie in the nonlocality of the expansion. It is important 
to understand as well how the singularities of the expan- 
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sion relate to symmetry breaking and the emergence of 
Goldstone modes which we will now discuss. 



V. ASYMPTOTIC FREEDOM 

The properties of branch cuts are clearly important in 
the partition function zeroes singularity picture as they 
relate to the emergence of a form of geometric bound on 
the modulus and argument of the zeroes in the complex 
expansion parameter plane |23l | . It is only the property of 
the zero nearest the real axis which is strictly important 
in determining the singularities in the N — >oo limit, but 
in the quantum case this also establishes the behaviour of 
an asymptotically defined source for the singularity. The 
scaling of the feature is therefore critically important to 
determining whether or not singularities of the zeroes in 
the spatial lattice polynomial expansion correspond to 
fixed points of this renormalization group equation de- 
fined along the real axis in (3. The sense in which we as- 
sociate the loop operator partition function zeroes treat- 
ment to the renormalization group equation approach is 
perhaps therefore closer to discussions of the strong in- 
teraction and a momenta dependent treatment Al- 
though the spatial volume is clearly important we are not 
treating the lattice spacing dependencies of this extent 
in a direct manner, and so the N — > oo properties of the 
zeroes are in a sense disconnected from the asymptotic 
treatment of the expansion parameter : in the sense that 
we would be interested to relate the correlation length 
and bare parameters of the system to finite N scaling 
properties of the volume. 

Never the less the symmetry breaking mechanisms of 
the system can be understood tentatively with relation 
to the Euclidean time extent and correlation length. If 
we take a conventional renormalization group analysis of 
a system like the 2-d 0(3) model we can understand that 
the correlation length, relating to the conventional renor- 
malization group /3-function coefficients, diverges expo- 
nentially in the limit of vanishing lattice spacing and so 
there are massless particles [J. However, for this par- 
ticular system, the Mermin-Wagncr-Coleman theorem 
prevents the massless Goldstone bosons interacting and 
so the model has a nonperturbatively generated mass 
gap [25| . Conversely other systems can be pictured in 
which the correlation length remains infinite and we have 
a broken symmetry, but reducing the Euclidean time ex- 
tent f3 eventually recovers a symmetric phase with a fi- 
nite correlation length. If the Euclidean time extent is 
very much smaller than the correlation length we can 
consider that all these systems undergo dimensional re- 
duction : we can treat one direction as effectively static. 
This then allows the continuum limit of these models to 
be defined with (3 — > 00. This limit is of interest be- 
cause we are interested in the asymptotic properties of 
[3 : which fixed point, or singularity on the real axis, we 
end up at in this limit. Notice that we have not aimed 
to fix at any point an association of the loop operator 



formalism with this dimensional reduction. In essence 
we therefore invert the above argument : having deter- 
mined singularities in the limit N — ► 00 we then consider 
the properties of these points on the real (3 axis where 
asymptotic convergence to the continuum is in principle 
guaranteed by construction. This is particularly useful 
form for the direct analysis of numerical lattice results 
since it is possible to understand from direct examina- 
tion of the singularities whether or not they correspond 
to genuine fixed points of a given model or are simply 
lattice spacing dependent cutoff effects. 

Considering the arguments defining a renormalization 
program for QED to all orders pcrturbativcly in terms 
of momenta, one must at some point make momentum 
cutoff dependent subtractions in the diagrams £|. Prc- 
scriptively constants are chosen so that the relevant Ward 
identities are satisfied at one particular value of momen- 
tum, and such that the renormalization functions do not 
contain singularities for the zero electron mass case. Sim- 
ilarly by construction the renormalized amplitudes of the 
diagrams do not contain infrared divergences for zero 
photon mass and nonzero electron mass. This prescrip- 
tion for infrared divergences of Gell-mann and Low's dis- 
cussed corresponds to the asymptotic limit of momentum 
cutoff being much bigger than the mass, and the solutions 
of the renormalization group equations for all momentum 
cutoff values and zero photon mass necessarily then de- 
fine a scale invariant field theory. One can break this 
scale invariancc in one of two ways ; by letting the mass 
be nonzero or by choosing a solution of the renormal- 
ization group equations at zero mass that is not a fixed 
point. One can therefore understand by analogy that for 
the loop operator approach all singularities in the spatial 
N extent can be connected to a finite (3 value with an 
asymptotic form, but that these points don't necessarily 
correspond to fixed point solutions if they do not have 
any support. To elucidate slightly : if a zero has an edge 
singularity associated with boundedness then its form is 
akin to this former case and we would have a broken 
symmetry with a correlation length that diverges as the 
continuum limit is approached, and the converse is that 
of a finite correlation length with a valid continuum limit. 
The third case is that of the mass gap, which would be 
otherwise under-defined in the analysis of say the case 
for example of the Ising model with external field, since 
some points of the complex phase plane diagram are nec- 
essarily undefined with the Laplace transform of the poles 
having essentially an under-defined contour |22| . 

Consider as an example the Haldane conjectured anti- 
ferromagnetic 0(3) quantum spin chains, where integer 
spins have a mass gap and half-integer spins are gap- 
less |l|- The low-energy effective theory for this system 
is the 2-d 0(3) model at vacuum angle 9 = tt, which is in 
the universality class of the k = 1 Wess-Zumino-Wittcn 
model. Numerical studies of this model indicate that the 
gap which is present for all 9 ^ tt, vanishes at the specific 
value 9 = tt. When the Euclidean time extent (3 is made 
finite this topological term also disappears. One can con- 
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elude from this system therefore that one would expect 
to see a singularity on the (3 axis from the loop opera- 
tor formalism, but that this should not scale directly to 
the continuum because the finite volume dependence of 
the scheme essentially introduces a divergent correlation 
length and one sees otherwise a pseudo-critical scaling. 
However the lattice singularity still exists on the finite 
volume and its relation to the asymptotic scaling can 
be measured directly from the lattice result, since again 
one can obviously study the gap dependence on finite sys- 
tem size directly in numerical simulation. The advantage 
with the new loop operator method is firstly that we can 
quantify the directly the asymptotic scaling properties of 
the gap itself, but also understand the evolution of other 
cutoff dependent quantum phase transition effects from 
the scaling arguments of related singularities. 

Preliminary results for the loop operator zeroes scheme 
are presented for the case of mixed integer - half-integer 
spin chains which share a related gapless property to the 
inte ger chains but with an additional saddle point ef- 
fect |2q . These admit a natural probabilistic loop oper- 
ator treatment as well following the pedagogical S=l/2 
AFM treatment schemes we have followed and discussed. 
An exact treatment has been recently given for the di- 
mensional reduction aspects of plaqucttc defined opera- 
tors over an extended phase space for the case of gauge 
field theories. This raises then raises the possibility of 
performing a similar analysis to the loop operator treat- 
ment for the case of gauge fields, as opposed to quasi- 
classical treatment, but the issues of numerically sim- 
ulation are by no means trivial to address. A general 
pathology of the factorised loop-cluster method, where 
the Monte Carlo decisions arc taken plaquette wise, is the 
statistical limitations of building up a long tailed proba- 
bility distribution, and extending to the case of gauge or 
complex valued fields one could then have also non posi- 
tive definite operator elements which do not yield a prob- 
abilistic interpretation at all. One option is to use a pro- 
jection operator basis acting on a conventional globally 
updated space to define pseudo-plaquette states [23 [28| . 
Strictly from the Laplace transform arguments we have 
reviewed and presented the only necessary criterion for a 
plaquette term on the space is that it contain at most one 
pole and have an continuous probability density function. 
In QCD, for example, a basis of related operators can be 
compiled over a given lattice gauge field ensemble such 
that the poles of the basis suppress unphysical state con- 
tributions through a natural ranking |29l | . Here we define 
first a Euclidean time dependent polynomial form for the 
global update scheme, independent of the related oper- 
ator basis, and introduce a limiting product such that 
the elements of these pseudo-plaquette operators contain 
at most one zero of the polynomial, which is easily done 
since the probability density functions of the operators 
on the globally generated gauge ensemble become uni- 
versally distributed at some critical scale [3(j. One can 
understand this latter result from general consideration 
of the central limit theorem in regard to defining the con- 



tour for the above Laplace transforms, and its relation to 
the stationary phase approximation in this regard |3lj . 

In lattice field theory calculations one can avoid some 
of the uncertainties associated with truncating a pcrtur- 
bative expansion in determining constants in the renor- 
malization group equations by defining various nonper- 
turbative schemes [32|[33. Constants are set through 
the chiral Ward identities for the lattice Green functions. 
For example with RI/MOM the conditions are set on 
Green functions at large external momenta, while in the 
Schrodingcr Functional scheme they are set in small lat- 
tice volumes. The lattice renormalization scales have to 
be larger than Kqcd, to connect with perturbation the- 
ory and smaller than the inverse lattice spacing to avoid 
cutoff effects, which is perhaps computationally expen- 
sive. With the RI/MOM approach as well the pseu- 
doscalar Green function, which couples to the Goldstonc 
boson, has a leading contribution diverging in the chiral 
limit. Similarly in building candidate supcrsymmetric 
lattice models one is interested in the Ward identities in 
the continuum limit in terms of renormalized perturba- 
tion theory |34j . Schemes involve, to finite order in per- 
turbation theory, adjusting a coefficient in local terms in 
the action but here as well it is possible for the action 
to become complex j3Ej. In both cases it is possible to 
envisage calculating a simple limiting operator factorisa- 
tion with some Euclidean time extent polynomial scheme 
defined such that the asymptotic properties of these lat- 
tice renormalization properties could be further studied 
directly in this context. 



VI. SUMMARY 

We have discussed, in the context of recent ideas on 
the role of dimensional reduction in quantum systems 
and reviews of loop operator methods a formulation of a 
partition function polynomial zeroes treatment for quan- 
tum systems. We have argued that in general the quan- 
tum case should be treated separately from the classi- 
cal case through a local projection operator argument 
for the factored product, suggesting the evolution of in- 
frared divergences is otherwise implicit to the quantum 
case. We have relaxed the condition that the expansion 
variable be expressed through an exponential function to 
this end, whilst still retaining its asymptotic structure. 
By working over an extended phase space we have then 
identified a limiting form with a connection to the con- 
tinuum, and the fixed point solutions of rcnormalisation 
group equations, using recent results for the treatment 
of exact probabilistic lattice dynamics. Related to this 
continuum limiting process we have then discussed the 
prospects for further evaluating gauge field theories for 
the purposes of perturbation theory independent renor- 
malization group analysis, and the treatment of theories 
with a non positive-definite source terms in the Euclidean 
action. The outline is necessarily general but further 
investigation of the limiting analyticity of the partition 
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function is proposed, with three cases with relevance to for future analysis and categorisation, 
symmetry breaking mechanisms envisaged as pertinent 
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